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Abstract 

Brown (1980, 1981) proved that the renewal function is concave if the inter-renewal dis- 
tribution is DFR (decreasing failure rate), and conjectured the converse. This note settles 
Brown's conjecture with a class of counter-examples. We also give a short proof of Shan- 
thikumar's (1988) result that the DFR property is closed under geometric compounding. 
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1 Introduction 

Structural relationships between the renewal function and the underlying distribution are of 
great interest in renewal theory. For a renewal process with decreasing failure rate (DFR) 
inter-renewal times, it is known that the renewal function is concave (Brown 1980). Conversely, 
Brown (1981) conjectures that DFR inter-renewal times are also necessary for the concavity 
of the renewal function. As shown by Shanthikumar (1988), there exist counter-examples to a 
discrete analogue of this conjecture. Brown's conjecture in the continuous case, however, has 
remained open. See Szekli (1986, 1990), Hansen and Frenk (1991), Shaked and Zhu (1992), 
Kijima (1992), and Kebir (1997) for related results and discussions. Also relevant is the work of 
Lund, Zhao and Kiessler (2006), who use hazard rates and renewal sequences to study reversible 
Markov chains. 

In this note we construct absolutely continuous distributions that do not have decreasing 
failure rates but nevertheless lead to concave renewal functions. That is, we give a definite 



answer to Brown's question in the negative. Our counter-examples have the following feature. 
On [0,ti] for some t\ > 0, the inter-renewal time has a decreasing failure rate; on [£1,^2] for 
some ti > t±, the failure rate strictly increases before decreasing again on ^2,00). It is shown 
that, for a suitable class of such distributions, if the increase in failure rate on [ti, ^2] is small 
enough, and the decrease shortly after t<i is fast enough, then the resulting renewal density is 
decreasing, i.e., the renewal function is concave. Section 2 presents the precise statements and 
illustrates with a numerical example. Section 3 contains the proofs. 

The renewal process is closely related to compound geometric random variables. In Section 4, 
by adapting the arguments of de Bruijn and Erdos (1953), we give an alternative proof of 
Shanthikumar's (1988) result that the DFR property is closed under geometric compounding. 

2 Concavity of the renewal function 

Let F(t) be a distribution function on R + = [0, 00) with F(0) = 0. Then the renewal func- 
tion M(t), i.e., the average number of renewals in [0,t], for a renewal process with underlying 
distribution F is given by 



(Some authors define M{t) + 1 as the renewal function; our results work with either definition.) 
If F(t) is absolutely continuous with density f(t), then so is M(t), and a version of its density, 
m(t), satisfies 



A positive function g(x), x € R+, is log-convex if logg(x) is convex on R + . A distribution on 
R + has DFR (decreasing failure rate), if its survival function is log-convex on R + . We recall 
two fundamental results relating M{t) to F{t). 

Theorem 1 (de Bruijn and Erdos (1953); Brown (1980); Hansen and Frenk (1991)). We have 

1. If F{t) has a log-convex density f(t), then the renewal density m(t) as in |I]) is also 
log-convex. 

2. If F(t) is DFR, then M(t) is concave. 

The question raised by Brown (1981) may be formulated as follows. 
Conjecture 1. If the renewal function M(t) is concave on R + , then F(t) is DFR. 





(1) 



2 



Shanthikumar (1988) resolves a discrete version of this conjecture by constructing a counter- 
example using auxiliary results on discrete Markov chains. It has also been noted that the 
discrete example does not generalize and the continuous case is still open. Our main result 
(Proposition [T]) finally disproves Conjectured) 

Proposition 1. Let < t\ < oo. Let f(t) be a density function that is positive on R + and 
continuously differ entiable on each of If., < k < 3, where Iq = [0, t\], L\ = [ti , £2] 5 I2 = 
[£2,^3], ^3 = [£3,00), with t 2 ,t 3 to be determined. That is, f(t) is continuous on R+, but f'(t) 
may jump at t^, k = 1,2,3. Assume the corresponding hazard rate function r(t) satisfies 

(i) r'{t) <0, te J ; 

(ii) on L\ we have 

X_ 

1 — ee 

for some e € (0, 1) where A > is determined by e and r{t\); 



r ( t ) = ~, ~v h<t<t 2 , (2) 



(Hi) on L 2 we have 

r'(t) < r\t) - /(0)r(t), t 2 < t < t 3] (3) 
(iv) r(t 3 ) < r(ti) and r'(t) < 0, t G I 3 . 

Then for small enough e > and t 2 — t\ > 0, both depending on the specification of r(t) for 
t € Lq only, the renewal density m{t) given by (CJ) decreases on R + . 

Note that r(t) strictly increases on [ii , £2] • Proposition Q] therefore settles Conjecture Q] in 
the negative. An example of a survival function F satisfying Conditions (i) and (ii) is 



Fit) 



i(e-* + l), 0<t<tx, 

( 4) 

ae~ xt - (3, ti <t< t 2 , 



where f3 > 0, and a and A are determined by j3 via 

F(t 1 +)=F(t 1 ) and F' (t 1 +) = F' (h-) . 

Specifically, 

A = [1 + (1 + 2/3)e tl ] ~\ a = (2A)~ V^ 1 )* 1 . 

The e in d2J corresponds to f3/a. Condition (iii) says that the hazard rate should decrease fast 
shortly after t 2 . An example based on (Jl]) that satisfies this condition is 

r(t) = r(t 2 )e^~^ 2 , t 2 < t < t 3 , 
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Figure 1: Illustration of a counter-example given by (JH), © and Q with t\ = 1, t 2 = 1.5, £3 = 2 
and = 0.02. 

which leads to 

F(t) = F{t 2 ) exp \-2r(t 2 ) (l - e^"*)/ 2 ) j , t 2 < t < t 3 . (5) 



For Condition (iv), we need £3 > t 2 + 2 log(r(i2)A'(£i)) to ensure r(t3) < r(ti), but r(t) can stay 
flat on t S ^3, which gives 

F(t) = J F(t 3 )e r(t3)(t3 -* ) , t>t 3 . (6) 

As an illustration, Figure 1 shows the survival function, density, hazard rate, and renewal 
function for a distribution as specified by (j3|), © and ([6]) with t\ = 1, t 2 = 1.5, £3 = 2, and 
/? = 0.02. The almost imperceptible decrease of m(t) on t E [1,1-5] is verified numerically as 
Proposition [1] guarantees monotonicity of m(t) for small enough (3 > and t 2 — t\ > but does 
not specify how small /3 or t 2 — t\ has to be. 

3 Proof of Proposition [I] 

We first establish a simple but useful identity. 
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Lemma 1. Let r(t) and F(t) denote the hazard rate and survival function, respectively, for a 
distribution with density f(t) on R + . Assume f(t) is absolutely continuous and f'(t) is bounded 
on every compact sub-interval o/R+. Then the renewal density m{t) as defined by HP satisfies 

m'(t) = r'(t)F(t) + / m'(x)[r(t - x) - r(t)]F(t - x)dx, t > 0. (7) 
Jo 

A discrete version of ([7]) can be traced back to Kaluza (1928). 

Proof of LemmaU\ The conditions guarantee that m(t) is absolutely continuous. In fact, we 
may differentiate under the integral sign in (pTJ) and get 

m '(t) = f'(t) +m(t)/(0) + I m(x)f'(t - x)dx, t > 0. (8) 

Jo 

Integration by parts then yields 

m'(t) = f'(t) + m(0)f(t)+f m'{x)f(t - x)dx, t > 0. (9) 

Jo 

We also have 

m'(x)F(t -x)dx = m(t) - f(0)F(t) - / m(x)f(t - x) dx 

Jo 

= f(t)-f(0)F(t), (10) 

where the first step uses integration by parts and the second uses (pQ). The identity ([7]) follows 
by expanding its right hand side and applying (flOl) and ([9]) to simplify. □ 

Proof of Proposition^ Since f'(t) is piecewise continuous, so is m'(t) as seen from ([8]). We have 
m'(0+) = r'(0+) < 0. Suppose m'(t) ever becomes nonnegative on Iq. Then letting t* be the 
smallest t £ (0, t±) such that m'(t) > we have m'{x) < 0, < x < £*, and by Condition (i) 
r'(t*) < 0, r(t* — x) — r(t*) > 0, < x < t*. It follows from ([7]) that m'(t*) < 0, a contradiction. 
Thus m'(£) < 0, i E Jo- I n fact, applying (J7J) again yields 

m'(t) < r'(t)F(t), t £ J , (11) 

where the left derivatives are used if t = t\. 
By ([2]) we have 

r'{t) = {\- l r{t) - l)r(t), *i < t < t 2j 

where A is determined from e via r (ti) = A/(l - ee Atl ). For fixed r{t x ) as e J. we have A f r{t x ) 
and hence r'(ti+) X 0. Denoting 5 = m'(t\—) — r'(t\—)F(t\), and noting 5 < by (fTT|) . we get 

m'(*i+) = r'(t 1 +)F(ti) + <5 < 0, 
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for small enough e > 0. Because m'(t) is continuous on (ti,^), and m'{ti+) < 0, we have 
m'it) < 0, t 6 (ti,i2), if £2 — h is small enough. Thus m(t) decreases on I\. 

Also, m(t) must strictly decrease on Ii- Assume the contrary and let t* be the smallest 
t € [t 2 ,t 3 ) such that m'(t+) > 0. Then © gives 

< m'(i*+) < f'(t*+) + m(0)/(t*), (12) 

because inside the integral m'(x) < 0, x € [0, £*). However, by ([3]) we have 

/'(**+) + m(0)/(t*) = F(t*) [r'(i*+) - r 2 (£*) + /(0)r(i*)] < 0, 

which contradicts (|12p . 

Finally, we show that m(t) decreases on I3 by applying ([7]) again. The assumptions rfe) < 
r(ti) and r'(t) < 0, t € ^3, ensure that r(t — x) > r(t), < x < t, t > t 3 , with strict inequality 
if t — x < t\. It is already shown that m'(x) < for x < £3. Thus ([7]) implies m'{t^+) < 0. The 
same argument proving m'(t) < for t E Iq then shows that m(t) decreases on I3. □ 

4 Preservation of DFR under geometric compounding 

Compound geometric random variables appear naturally in areas such as queuing theory (see, 
e.g., Szekli 1986) and financial risk modeling. It is well known that log-convexity is closed under 
geometric compounding (this is essentially Part 1 of Theorem [T]). Shanthikumar (1988) showed 
that the DFR property is also closed under geometric compounding. This was achieved by 
establishing auxiliary results on discrete Markov chains. It may be worthwhile to note that the 
argument of de Bruijn and Erdos (1953) can be adapted to give a short proof of Shanthikumar's 
(1988) result (Part 1 of Theorem [2]). The same argument yields a parallel result (Part 2 of 
Theorem [2]) concerning the increasing failure rate (IFR) property. 

Theorem 2. Let X be a random variable on N = {1,2,...} and let T be a geometric with 
parameter p € (0,1), i.e., Pr(T = n) = pq™' 1 , n = l,2, q = 1 — p. Define the random sum 
Y = Ylk=iXk where are independent (and also independent ofT) and identically distributed 
as X . 

1. 7/logPr(A > n) is convex in n € N, i.e., X is discrete DFR, then so is Y. 

2. If logPr(y > n) is concave in n £ N, i.e., Y is discrete IFR, then so is X. 
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Proof. Denote 

f n = Pv(X = n), F n = Pr(X>n), g n = Pr(Y = n), G n = Pr(Y>n). 
We have the recursions 

n— 1 n— 1 

9n=pfn + q^2fk9n-k, G n = F n + fkG n -k, n=l,2, .... (13) 

k=l k=l 

The following identity is analogous to Equation (7) of de Bruijn and Erdos (1953); Hansen 
(1988) uses similar identities for compound Poisson probabilities (see Yu 2009 for related work). 
It is proved by expanding the right hand side and then applying (113H . 

F n (G n+ 2G n — G^+i) =pG n (F n+ 2F n — F^ +1 ) 

n 

+ <?y^(-^n+l/fc-l - Fnfk){G n +lG n+ i-k — G n G n+2 -k)- (14) 
k=2 

In particular, G3G1 — G\ = p{F^,F\ — F 2 2 ). Assuming F n is log-convex, we get G n+ \G n+ i_k > 
G n G n+ 2-k, 2 < k < n, and G n+ 2G n > G^ +1 , n > 1, by induction from (fT4"|) . Thus G n is 
log-convex in n € N, i.e., Y is discrete DFR, and Part 1 is proved. Similarly, assuming Gn is 
log-concave, we get F n+ i/ fc _i < F n f k , 2 < k < n, and F n+2 F n < F% +1 , n > 1, by induction. 
Thus -F n is log-concave in n 6 N, i.e., X is discrete IFR, and Part 2 is proved. □ 
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